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Abstract. Context. Since its launch in November 2004, Swift has revolutionised our understanding of gamma- 
ray bursts. The X-ray telescope (XRT), one of the three instruments on board Swift, has played a key role in 
providing essential positions, timing, and spectroscopy of more than 300 GRB afterglows to date. Although Swift 
was designed to observe GRB afterglows with power-law spectra. Swift is spending an increasing fraction of its 
time observing more traditional X-ray sources, which have more complex spectra. 

Aims. The aim of this paper is a detailed description of the CCD response model used to compute the XRT 
RMFs (redistribution matrix files), the changes implemented to it based on measurements of celestial and on- 
board calibration sources, and current caveats in the RMFs for the spectral analysis of XRT data. 
Results. We show that the XRT spectral response calibration was complicated by various energy offsets in photon 
counting (PC) and windowed timing (WT) modes related to the way the CCD is operated in orbit (variation 
in temperature during observations, contamination by optical light from the sunlit Earth and increase in charge 
transfer inefficiency). We describe how these effects can be corrected for in the ground processing software. We 
show that the low-energy response, the redistribution in spectra of absorbed sources, and the modelling of the 
line profile have been significantly improved since launch by introducing empirical corrections in our code when it 
was not possible to use a physical description. We note that the increase in CTI became noticeable in June 2006 
{i.e. 14 months after launch), but the evidence of a more serious degradation in spectroscopic performance (line 
broadening and change in the low-energy response) due to large charge traps {i.e. faults in the Si crystal) became 
more significant after March 2007. We describe efforts to handle such changes in the spectral response. Finally, 
we show that the commanded increase in the substrate voltage from to 6 V on 2007 August 30 reduced the dark 
current, enabling the collection of useful science data at higher CCD temperature (up to -50° C). We also briefly 
describe the plan to recalibrate the XRT response files at this new voltage. 

Conclusions. We show that the XRT spectral response is described well by the public response files for line and 
continuum spectra in the 0.3-10 keV band in both PC and WT modes. 

Key words, gamma-ray: bursts - X-rays: general - Instrumentation: detectors - Methods: numerical - Mission: 

Swift 



1. Introduction 

Successfully launched on 2004 November 20, the Swift 
gamma-ray burst satellite (Gehrels et al. 2004) consists 
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of three instruments: the wide-field of view, gamma-ray 
burst alert telescope (BAT; Barthelmy et al. 2005) and 
two narrow field instruments (NFIs), the X-ray telescope 
(XRT; Burrows et al. 2005) and the UV/optical telescope 
(UVOT; Roming et al. 2005). Thanks to the unique abil- 
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ity of Swift to slew automatically after the BAT trigger, 
bursts are typically within the field of view of the nar- 
row field instruments within a couple of minutes after the 
trigger. Thus, the XRT routinely provides positions with 
an accuracy of a few arc-seconds to the GRB community 
world-wide and measures the early X-ray light-curves and 
spectra of most GRB afterglows at which it is promptly 
pointed. Up to 2008 August 11, a total of 357 GRBs were 
detected by the BAT, of which 314 GRBs were observed 
by the XRT, although only 297 X-ray afterglows were thus 
detected by the XRT. Of these burst observations, 253 
were prompt slews (i.e. less than 300 s after the BAT trig- 
ger), and the XRT detected 94%[]of those bursts (i.e. 238 
afterglows) . 

The Swift-XBT spectro-temporal observations associ- 
ated with on-board and ground-based multi-wavelength 
observations have actively driven GRB science in the past 
three years by shedding new light on the physics of these 
objects. X-ray observations have revealed previously un- 
expected behaviour, including: multiple temporal breaks 
observed in some GRB light-curves inconsistent with the 
standard afterglow models (e.g. Zhang et al. 2006, Nousek 
et al. 2006, O'Brien et al. 2006, Willingale et al. 2007) 
and the discovery of X-ray flares observed in about 50% 
of the Swift afterglows (e.g. Falcone et al. 2006, 2007; 
Chincarini et al. 2007; Goad et al. 2007; Godet et al. 
2007a, 2006; Burrows et al. 2007). Recent work has shown 
that some spectra cannot be fit by simple absorbed power- 
laws due to curvature, which can be interpreted as spectral 
breaks caused by the temporal shift of the energy peak 
through the XRT energy band during X-ray flares (e.g. 
GRB 051117A: Goad et al. 2007 and GRB 050822: Godet 
et al. 2007a). In some other cases, the presence of an extra 
component was suggested (see Butler 2007, Moretti et al. 
2008a for a general study) such as a blackbody compo- 
nent, which could be interpreted as: possible emission of 
a jet cocoon; the first-ever detection of a shock breakout 
from a massive star, likely a Wolf-Rayet, in the peculiar 
event GRB 060218 (Campana et al. 2006); possible pho- 
tospheric emission for some X-ray flares in GRB 050822 
(Godet et al. 2007a). 

More detailed X-ray spectral analyses have shed light 
on the environment of GRBs and the nature of the pro- 
genitor itself. Campana et al. (2007) showed evidence that 
the progenitor of the high redshift GRB 050904 at z = 6.3 
was located in a dense molecular cloud with a metallic- 
ity Z > 0.03 Zq. Using detailed modelling of the low- 
energy part of the X-ray spectrum, Campana et al. (2008a) 
showed that the progenitor of GRB 060218 was likely to 
be a massive star characterized by a fast stellar rotation 
and initial sub-solar metallicity giving, for the flrst time, 
direct evidence about the properties of GRB progenitors. 

Thanks to its large energy band-pass and the abil- 
ity to rapidly schedule targets. Swift is also a power- 
ful tool for obtaining essential insights into the physics 



of non-GRB objects, especially transients. For example, 
the XRT has provided essential spectroscopic information 
since the start of the outburst of the recurrent nova RS 
Ophiuchi in 2006, enabling, for the flrst time, detailed ob- 
servations of the evolution of the Super- Soft Component 
in such an object (e.g. Bode et al. 2006, Hachisu et al. 
2007). It also provided important spectroscopic informa- 
tion during the follow-up of several transient events (e.g., 
blazars: Tramacere et al. 2007; X-ray binaries: Esposito et 
al. 2007, Wijnands et al. 2007, Rykoff et al. 2007, Romano 
et al. 2007, Brocksopp et al. 2005; Comet 9P/Tempel 1: 
Willingale et al. 2006; active galactic nuclei: Molina et al. 
2007, Tueller et al. 2007). From the beginning of 2008, 
the fraction of time spent on non-GRB science targets has 
been ~ 45%, and this will increase in the future. 




^ see http://swift.gsfc.nasa.gov/docs/swift/archive/grb-table/ 
grb- stats, php 
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Fig. 1. Energy resolution (a = ^^™ ) versus energy as 
measured during the ground calibration at the Leicester 
calibration facility using the in-flight CCD-22 cooled at 
-100° C. The resolution was measured over the central 
200 X 200 pixel window. The curve corresponds to the 

theoretical prediction given by a = uj^ ^ -\-r^ where F, 
uo and r are the Fano factor, the photon energy, the 
pair creation energy coefflcient and the readout electronic 
noise, respectively. The thermal noise could be neglected 
with respect to the readout electronic noise due to the low 
temperature used during the measurements. 

The results above show that X-ray spectroscopy is a 
powerful and essential tool for obtaining physical insights 
for celestial objects. These results depend on the intrin- 
sic performance of the detector and the calibration of its 
spectral response. The XRT uses a front-illuminated e2v 
CCD-22, offering good spectroscopic performance with an 
energy resolution (FWHM) before launch of 135 eV at 
5.9 keV (see also Fig. [l]). The CCD-22 is a three phase 
frame transfer device, which utilises high resistivity (typ- 
ically 4000 Vt cm) silicon and an open electrode structure, 
originally designed for, and used in, the EPIC (European 
Photon Imaging Camera) MOS cameras on-board XMM- 
J Newton. The open electrode structure signiflcantly im- 
proves the quantum efflciency of the device at low energy. 



Godet et aL: Modelling of the Swift-XRT spectral response 



3 



N,, = 1,2400 



Energy Channel CH 
Channel width = 5eV 



N = L N„ 



X-ray energy E 
randomised in CH 



Random position in a pixel 
40^ X 40^m (X,Y) 



Check if the photon passes through the 
electrode, finger or buried channel structure 



Input parameters 

Linear Abs. coeff. 
(Si, SiO^, Si^N^) 

Nj^j.: nmn of photons 
per channel 

CTE_: Readout noise 



V 



Absorption a 


t depth z using 


Prob(z) = 


exp(-z/M,(E)) 



Si Ktt fluorescence 



(Variance on Qg 



Electron cloud generation 
r^ = 0.0044 E^i ^5 fmi 

Qn = E/w 



E>Ei^ = 1.839 keV 
and 
Proba < 0.044 



1 ^ 

X 



= E-1.74 keV 





I 


Depletion deptli 

I'd 





z >Z^ 



TRUE 



/tjeneration of*^ 
_ ' a Si Ka photon * 
I withX%Y%z' J 
\e, = 1.74 keV/ 



FALSE 



Field-free region 



JrRUE| 



Substrate region 

r. 



Collected charge 



T 



Charge entering 
depletion depth 

% 



Charge entering 
Field-free region 



Surface charge loss 
Qtd=Q(rt)*F(^.E) 
F(z,E) < 1 



5x5 Matrix (PC) 
5 X 1 Matrix (WT) 
Q(n.m) 



Readout noise I 



i 


Pulse Height distributions 
normalised by N^^^ 
for channel CH 


Q(n,m) > Threshold 
Event grade selection 







Fig. 2. Architecture of the CCD response model used to compute the XRT RMFs. The RMFs contain 2400 input 
energy channels with a 5 eV energy width. In each energy channel, N^ax > 5 x 10^ photons with an energy randomised 
in a given bin are considered in order to avoid adding numerical noise to the computation. jiz^E) corresponds to 
the linear absorption coefficient in the silicon, which is a function of the energy of the incoming photons {E)^ while 
CTEg^p corresponds to the charge transfer efficiency during the readout process in the serial and parallel direction, 
respectively (see Section 3.2.2). Zd, zjj correspond to the depth of interaction of the incoming photons {z = 
corresponding to the surface of the Si active bulk), the depletion depth (an input parameter) and the field- free region 
depth defined as zjj = Azccd — — (where Zg = 200 jam is the thickness of the substrate and Azccd = 280 /im is 
the thickness of the CCD). F = 0.115 and co = 3.68 eV per pair electron/hole are the Fano factor and the pair creation 
energy coefficient in the silicon, respectively, while Qo corresponds to the initial charge generated by the interaction 
of an incoming photons within the device. The different equations giving the total charge spreading radius (rt) in the 
different layers of the CCD and the collected charge Q(rt) can be found in the Section 2.1.1| f{E^z) corresponds to a 
function describing the charge losses at the surface of the CCD (see Section 2.1.3). It depends on energy and depth of 
interaction of the incoming photons. Mostly photons interacting in the opened electrode area will suffer from surface 
charge losses (see the text in Section2.1.3). The dashed lines correspond to the case where a 1.74 keV Ka fluorescence 
photon is emitted following the interaction of an incoming photon with an energy E > E^^ = 1.839 keV where E}y is 
the electron binding energy for the Si K shell. The probability of generating such a fluorescence photon is 4.4%. 



The CCD imaging area consists of a 600 x 600 pixel array 
with a pixel size of 40 x 40 jim^ . Forty microns corresponds 
to 2.36 arc-seconds in the XRT field of view (e.g. Short et 
al. 2002). To mitigate the effects of pile- up, the CCD can 



automatically switch between different readout modes de- 
pending on the source brightness once the spacecraft is 
settled on the source (Hih et al. 2004, 2005): 
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1) Photo-diode (PD) mode at the highest count rates 
with a 0.14 ms time resolution and no spatial infor- 
mation; 

2) Windowed Timing (WT) mode at moderate count 
rates, which uses a 200 column window covering the 
central 8 arcminutes of the XRT field of view and pro- 
vides 1-D spatial imaging information. In that mode, 
the columns are clocked continuously to provide tim- 
ing information in the trailed image along each column 
with a 1.8 ms time resolution, at the expense of imag- 
ing information in this dimension (pixels are binned by 
a factor of 10 along columns); 

3) at lower count rates. Photon Counting mode (PC), 
which provides 2-D spatial imaging information, but 
with a 2.5 s time resolution. 

The PD mode has not been used since the CCD was hit by 
a particle (micro-meteoroid) on 2005 May 27 (Carpenter 
et al. 2006), due to the apparition of several bad columns 
on the detector, so we will focus only on the PC and WT 
modes in the remainder of the paper. 

The XRT effective area (EA) is made up of the re- 
sponse of the mirrors, the filter placed in front of the CCD 
detector to reduce optical and UV photon loading and the 
CCD detector response. The mirror and filter responses 
are included in the auxiliary response files (ARFs), while 
the CCD response is included in the redistribution matrix 
files (RMFs). 

To avoid any misinterpretation in the spectral anal- 
ysis, it is essential to have a good understanding of the 
instrument response and its operational limitations. The 
aim of this paper is to describe in detail our CCD response 
model used to compute the RMFs, the changes that have 
been made to the RMFs and ARFs since the launch (based 
on in-flight calibration) and the caveats to be aware of in 
the spectral analysis of XRT data when using the cur- 
rent (version 11) response files distributed in the CALDB 
release on 2008-06-25. An overview of the other aspects 
of the XRT calibration can be found in Campana et al. 
(2008b) and a detailed study of the XRT background can 
be found in Moretti et al. (2008b). This paper is organised 
as follows: 

2 - We describe the physics implemented in our RMF 

code. We also briefly describe the spectroscopic per- 
formance of the RMFs in both PC and WT modes 
prior to launch using ground calibration data. 

3 - We discuss the in-flight operation of the XRT and its 

impact on the calibration of the XRT response, the 
calibration program, and in detail post-launch changes 
made to the CCD spectral model (low-energy response, 
line shoulder, shelf, RMF redistribution) and the ARFs 
based on celestial target calibration. 

4 - We present the spectroscopic performance of the XRT 

for several celestial targets compared to observations 
with other X-ray instruments, as well as caveats (line 
broadening due to the build-up of charge traps on the 
CCD and changes in the in-orbit operation of the CCD 
caused by raising the substrate voltage from V to 6 V 



in 2007 August 30) for the spectral analysis when using 
the current RMFs and ARFs. 
§5 - We present the main conclusions of the paper. 

2. Computation of the spectral response 

The XRT spectral response is based on a physical model 
describing the interaction of photons within the silicon 
bulk of a CCD-22 pixel via Monte-Carlo simulations. A 
Monte- Carlo simulator was developed at the University of 
Leicester in order to generate the RMFs for all the XRT 
readout modes (Short et al. 2002, Mukerjee et al. 2004). 
The CCD response model was refined by implementing 
ad-hoc corrections in our code based on the ground and 
in-flight calibration. We discuss here the main steps used 
to generate the RMFs. 

2.1. RMF generation 

Figure [2] shows the architecture of our CCD response 
model. To compute the RMFs, we stack simulated spec- 
tra of monochromatic X-rays using simple geometry of 
a CCD-22 pixel as input parameters, including the elec- 
trode structure, the finger and the open area (see Fig. [3| 
and the linear absorption coefficients of the different CCD 
compounds. The linear absorption coefficients used in our 
simulations include the effects of the rapid oscillatory X- 
ray Absorption Fine Structures (XAFS) just above the 
absorption edges of the CCD constituents (mainly N, O 
and Si) as measured by Owens et al. (1996a) for the JET-X 
CCD (see Fig[4|. It is essential to take them into account 
to obtain proper modelling of the spectral response around 
the edges. 

2.1.1. Carrier generation and charge spreading 

0.2-10 keV photons impinging on the detector have a high 
probability of interacting via the photoelectric effect. This 
results in the generation of a charge cloud with an initial 
charge Qo given by: 

Qo = e X + uo X (1) 

where F and uj correspond to the Fano factor and the pair 
creation energy coefficient in the silicon, respectively. R is 
a random number uniformly distributed between -1 and 
1; e is the elementary charge of the electrons (e = 1.6 x 
-|^Q-i9 j£ ^Yie incoming photons have an energy E above 
El) = 1.839 keV (i.e. the binding energy of the K-shell 
electrons in the silicon), then there is a 4.4% probability 
of producing a Si Ka fluorescence photon at 1.74 keV. 
In this case, the residual energy \^ E^ = E — 1.74 keV. 
Otherwise, we have: E^ = E. 

The newly formed charge cloud has an initial 1 a radius 

of 

= 0.0044 El'^^ fim (2) 
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T3 = T2 exp(-0.25//gj ) exp(-0.3//gio^) 
nna = T^exp(-O.Vg.^^ ) exp{-0.4^gjQp 

^ = T3exp(-0.25p'gj ) exp{-D.4/rgjQ^) 



Fig. 3. Pictures of the CCD-22 structure and geometry of a CCD-22 pixel as used in our spectral response code. (Top 
left and right panels) Top and side view of a CCD-22 pixel, respectively. These images were obtained by a scanning 
electron microscope provided by e2v. The different parts of the pixel are shown on these two images. (Bottom left and 
right panels) Geometry of the CCD-22 pixel as used in our CCD response model (the dimensions are given in units of 
/im): (left) Top view of the simplified geometry of a CCD-22 pixel. The thicknesses of the compounds (Si, Si02 and 
Si3N4) present in the finger and electrode structures are shown in the equations giving the transparency of each part 
of the electrode structure on top of the silicon bulk (the coefficients T^=i^4); (right) Side view showing the different 
layers formed in the CCD (see text for more details). The epitaxial layer is 80 fim thick (p-type) and the substrate is 
200 /im thick. The in the z coordinate corresponds to the top of the Si active bulk, the electrode being on top of it. 
The incident X-ray photons enter from the top. 



assuming that the spatial distribution of the charge follows 
a Gaussian (Fitting et al. 1977). This charge cloud is col- 
lected in the buried channel in the depletion region after 
spreading in the bulk of the detector. The buried channel 
is defined as the potential well where the photo-generated 
charge will be collected and stored in a given CCD pixel 
between each readout. For the CCD-22, the buried channel 
covers the entire pixel size. The charge cloud may spread 
into more than one pixel, depending on its position with 
respect to the electrode structure and its depth, which is 
a function of the value of the linear absorption coefficients 
at the energy of the incoming photons (E) as follows: 

^ = --^ln(P + e) (3) 

where 2:, jJ^ziE) and P are the depth of interaction of the 
incoming photons, the linear absorption coefficient and a 



number uniformly and randomly distributed between 
and 1, respectively. A small offset e <C 1 is introduced to 
avoid any singularity when P = 0. 

The CCD pixel depth is divided into three regions: 

- The depletion region where the charges experience 
the full effect of the electric field (E). The radius of the 
charge cloud is modified as it drifts to the buried channel. 
The contribution to the charge cloud radius due to the 
spreading in the depletion region is given by: 

u = m ( (4) 

where Zd is the thickness of the depletion depth, Na - the 
silicon doping concentration {Na = 8 x 10^^ cm~^), Cg - 
the Si permittivity (e^ = 1.044 x 10"^^ F cm-^), k - the 
Boltzmann constant and T - the temperature. We use a 
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teraperature of — 60°C in our code, which corresponds to 
the average CCD temperature in-orbit (see Section 3.2.3). 
CCD temperature fluctuations by ±10°C do not change 
dramatically the r^^-values computed using Eq. |4] (i.e. by 
less than 2%). Equation. [4] is undefined when z = z^- 
However, this singularity is not physical, so we introduced 
a small offset e = 0.002 <C 1 to avoid this. All the charge 
reaching the depletion region is collected in the buried 
channel. 



100.000 
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0.100 
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Fig. 4. Linear absorption coefficients of the CCD-22 com- 
pounds: (solid line) Si; (dotted line) Si02; (dashed line) 
Si3N4. The small inset shows the XAFS just above the Si 
edge. The energy of the Si edge is slightly different depend- 
ing on the compound considered. This is due to differences 
in the crystal lattice structure. 



- The field-free region where the charge no longer feels 
the effect of the electric field (i.e. a region where E ~ 0). 
In this region, the charge experiences diffusion, recombina- 
tion and trapping in the device lattice, so that the charge 
cloud can spread into several adjacent pixels. The contri- 
bution to the charge cloud radius, due to spreading in the 
field-free region, is given by (Janesick et al. 1985): 




(5) 



where zjj is the thickness of the field-free region and 
Z2 = z — Zd > is the depth of the interaction in the 
field-free region. Assuming that the recombination losses 
can be neglected, the charge reaching the depletion layer 
boundary is given by (Hopkinson 1983): 



Qff = Qo- 



cosh 



/ Zff-iz-Zd) \ 



cosh 



(6) 



where Ljj is the value of the diffusion coefficient in the 
field-free region {Ljj ~ 4000 /im). Equations, [s] and [6] ini- 
tially used to describe the diffusion in the field-free region 
have been modified to follow the formalism described in 



Section 3 in Pavlov & Nousek (1999), because the assump- 
tion that the radial distribution of the charge is Gaussian 
is no longer valid in this region (e.g. Pavlov & Nousek 
1999; see also Section 3.3.3). 

- The substrate region in which most of the charge pro- 
duced is lost by diffusion and recombination, because the 
electron diffusion length is much shorter than in the two 
other regions. The contribution to the charge cloud ra- 
dius due to spreading in the substrate region is given by 
(McCarthy et al. 1995) 



(7) 



where Zs is the thickness of the substrate, Ls - the diffu- 
sion length in the substrate which is much shorter than 
in the field- free region (Lg = 10 /im for Na ~ 10^^ cm~^) 
and Zs = z — Zd — Zff is the depth of the interaction in 
the substrate. The diffusion length is replaced by Zs when 
using Eq. [7|to avoid having {zg — zs) > Lg (Short et al. 
2002). 

The total 1 a radius of the charge cloud reaching the 
buried channel is then given by the quadratic summation 
of the relevant radii assuming that the charge cloud profile 
due to radial diffusion is normally distributed: 




+ + + 



'ff 



(8) 



The contribution of can be neglected with respect 
to that from the spreading radii r^^, rjj and r^. 

Figure. [3] gives the thicknesses of these different lay- 
ers as used in the CCD response model. The thickness of 
the depletion depth is an input parameter in our model. 
From quantum efficiency (QE) measurements performed 
at the Leicester calibration facility, we demonstrated that 
a 27/im-thick depletion depth matches the on-ground cal- 
ibration data using a substrate voltage of Vgs = OV, as 
initially set in flight. 

2.1.2. Event recognition process 

The spreading of the charge cloud over several pixels im- 
plies that it is essential to use an event grade recognition 
process to retrieve the right energy of the incoming pho- 
tons. We use the same event grade recognition process 
as implemented in the ground analysis software (xrtdas 
software package]^ Capalbi et al. 2005). 

An event is formed when a pixel has an analogue-to- 
digital converted (ADC) charge greater than the event 
threshold (80DNr]). This threshold is used to avoid pix- 
els containing only noise from swamping the telemetry. 
In PC mode, the event data are telemetered as a 3 x 3 
pixel matrix centred on the pixel with the most charge. 
Only the surrounding pixels for which the charge is above 
the split event threshold (40 DN) are considered in the 
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Fig. 5. List of the XRT event grades: (left) for PC mode, a 3 x 3 pixel matrix centred on the pixel containing the 
most charge is used to define the event grade, while for WT mode (right) a 7 pixel vector is used. The grades 0-12 in 
PC mode and the grades 0-2 in WT mode are considered as valid events due to X-rays as shown by the grey boxes in 
the Figure. 



grade mapping, the other pixels being discarded. The split 
event threshold was fixed at 40 DN to minimise the in- 
crease of the sub-threshold charge losses which result in a 
wing on the low-energy side of the line profile, and hence 
to minimise the degradation of the energy resolution (see 
Sections 2.1.3 and 3.3.3). 



Since the WT mode only offers 1-D spatial imaging 
information (see Section 1), it is impossible to use a 3 x 3 
pixel matrix in WT mode to classify the X-rays events. 
Instead, a seven pixel vector centred on the pixel with the 
most charge is used. In WT mode, the central and split 
event thresholds are both set at 80 DN. Figure [5] shows 
the list of the 32 and 15 event grades which can be formed 
in PC and WT modes, respectively. The PC grades are 
similar to those defined for the XMM-Newton MOS data 
(Burrows et al. 2005). 



To eliminate events due to charged particles, and to 
obtain good energy resolution, we consider grades 0-12 
in PC mode and grades 0-2 in WT mode as valid X-ray 
events. From Fig. [5) the WT grade includes the PC 
grades 0, 1 and 3 (i.e. the mono-pixel and vertical split 
events in PC modes), while the WT grades 0-2 includes 
the good PC grades 0-12 as well as possible higher PC 
grades 15, 17, 19, 21, 23, 25-29. 

The pre-launch RMFs (v006) were released for three 
grade selections (0, 0-4 and 0-12) in PC mode and two 
grade selections (0 and 0-2) in WT mode. This was to 
offer the user the choice of higher spectral resolution at 
the cost of lower effective area. The decision was made 
after launch not to upgrade the PC grade 0-4 RMFs after 
the release of the v007 response files, since the PC grade 0- 
12 RMF offers a higher quantum efficiency at high energy 
and its calibration was sufficiently good. 
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2.1.3. X-ray spectrum from monochromatic radiation 

The X-ray event energy spectrum resulting from 
monochromatic radiation with an energy E significantly 
differs from a simple Gaussian. It consists of multiple 
components: a Gaussian photo-peak with a shoulder on 
the low energy side, a shelf extending to low energies, 
and, at the very lowest energies, the high-energy side of 
a noise peak above the on-board central event threshold 



(see Section 2.1.2). For photon energies above the Si K- 



shell edge (1.839 keV in the silicon bulk) two additional 
features appear: an escape peak of energy E — Esi, and a 
Si Ka fluorescence peak at Esi = 1.74 keV. 

The DN value of the minimum energy threshold 
{thres) in the spectra depends on the mode, the grade 
selection and the settings for the event threshold and the 
split event threshold in the ground processing code. The 
default threshold configuration for PC mode is 80 DN and 
40 DN, respectively, and for WT mode both are set to 
80 DN. Therefore, in PC mode, the default minimum en- 
ergy threshold for grade events is 80 DN while it is set 
at (80 -h 40) DN for double split events, (80 + (2 x 40)) DN 
for triple split events and (80 + (3 x 40)) DN for quadruple 
split events. In WT mode, the minimum energy threshold 
is set at 80 DN and 80 + 80 DN for grade and grades 1 
& 2 events, respectively. The t/ires-values once expressed 
in units of eV [i.e. thres x Cq x G where Cq and G (= 10) 
are the multiplicative DN to PHA (pulse height ampli- 
tude) gain factor and the global PHA to PI (pulse in- 
variant) gain factor, respectively - see Eq. 10 slightly in- 



creased over time, since the gain coefficient Cq has in- 
creased by about 3% from launch to June 2007 (see Fig. 12 



and Section 3.2.2). 



Mil Ka 
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Fig. 6. Spectra from the Mn Ka and KP lines (5.9 and 
6.5 keV) from an on-board ^^Fe radioactive source in 
WT mode (black crosses - grades 0-2) and PC mode (red 
crosses - grades 0-12). 

At low energy, a large fraction of the photons interact 
at the interface between the Si02 layer placed on top of 



the silicon bulk and the active silicon bulk itself. Short et 
al. (2002) measured the lost fraction of the total charge 
formed at that interface. This energy loss results in a low- 
energy wing in the line profile as well as a shelf extending 
to low energy. The effects of the surface charge losses start 
to become very significant below 0.5 keV, since the photo- 
peak disappears almost completely and the resulting line 
profile is shifted to lower energies (see Short et al. 2002). 
The authors proposed that the surface loss effect may be 
due to the charge-state of the surface oxide so that the 
oxide layer is charged enough in the etched regions (open 
electrode areas shown in white in the bottom left panel in 
Fig. [3| to cause a local turn-over of the potential (see Fig. 
7 in Short et al. 2002). In this case, the charges formed 
near the oxide layer will move to the surface rather than to 
the buried channel. The surface charge losses are a func- 
tion of the depth of interaction and the energy. To model 
this double dependency, the XRT energy range is divided 
into 12 energy bands {AEi=i^i2). In each energy band, 
the surface charge losses were initially defined by a set 
of linear functions /(A£^^=i^i2, z) depending on the depth 
of interaction z inside the silicon bulk, the coefficients of 
the function f{AEi=i^i2, z) being empirically derived from 
spectroscopic measurements performed at different ener- 
gies at the Leicester calibration facility using radioactive 
elements. 

At higher energies, when the photons interact more 
deeply in the CCD (above 2 keV), other processes start to 
become more important in the production of the shoulder 



and the shelf: (i) sub-threshold losses (see Section 3.3.3); 
(ii) diffusion, recombination and trapping in the bulk of 
the detector; (iii) inhomogeneity of the electric field in the 
depletion depth; these act in addition to the surface losses. 
The exact shape of the shoulder and the shelf depend on 
the readout mode (see Fig. [g] and the top panel in Fig.[7|). 

2.2. Performance of the pre-launch RMFs 

The XRT CCD was calibrated using the EPIC calibration 
facility at the University of Leicester, using sixteen dis- 
crete energies covering the 0.3-10 keV energy range. These 
measurements were used to interpolate and generate the 
redistribution component of the XRT RMFs. The RMFs 
were initially calibrated for a value of the substrate volt- 
age of Vss = V, which was the value used from launch 
to 2007 August 30 (see Section |43|. 

The pre-launch RMFs were computed to model the 
spectral response in the central region of the CCD, i.e. 
a window of 200 x 200 pixels (a field of view of about 
7.9 X 7.9 arc-minutes^). This area also corresponds to the 
area on the CCD where most of the GRB X-ray afterglows 
are located after the spacecraft slews. 

Although the pre-launch RMFs showed good agree- 
ment between the calibration data and the model in many 
circumstances, as displayed in the top panel of Fig. [T) there 
were still some modelling issues. The low-energy shoulder 
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Fig. 8. CCD-22 simplified schematic diagram. The four 
^^Fe corner sources permanently illuminating areas of the 
CCD outside the field of view are shown as grey circles. 
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Fig. 7. (Top panel) PC grade spectrum measured at the 
Leicester calibration facility from the Si Ka and Si Kp 
hues (1.74 and 1.84 keV) with background O Ka (0.525 
keV) and bremsstrahlung continuum, together with pre- 
launch RMF v006. (Bottom panel) PC grade 0-12 spec- 
trum measured at the Leicester calibration facility Mn Ka 
and KP lines (5.89 and 6.5 keV) from an on-board ^^Fe 
source together with the pre-launch PC RMF v006. The 
low-energy shoulder and shelf are not fully reproduced by 
the model. 



Many of our calibration targets are also observed 
by other X-ray observatories such as XMM-Newton^ 
Chandra, Suzaku and RXTE, enabling us to perform cross- 
calibration campaigns. So far, five cross-calibration cam- 
paigns with the XMM-Newton EPIC cameras on variable 
sources have been made (see Table [l] and Section 4.1 ). We 
also use the set of stable calibration sources observed by 
other X-ray instruments to compare and improve the per- 
formance of our spectral response (see Table [l] and Section 
4.1). 



of the high-energy {E > 2 — 3 keV) line profiles and the 
shelf needed improvement (see the bottom panel in Fig.[7|). 

3. In-flight calibration 

3.1. Calibration overview 

Regular in-flight spectroscopic calibration observations of 
a set of well-known celestial objects are performed (with a 
frequency of about 6 months) in order to monitor changes 
in the response (see Tablejl]). The fraction of time spent on 
calibration since the launch is ~ 7% of the total in-orbit 
time. In addition, we make use of four ^^Fe calibration 
sources which are located in each corner of the CCD and 
which permanently illuminate a small fraction of the CCD 
area outside the XRT field of view as shown in Fig. [Sj 
Before the focal plane camera assembly (FPCA) door was 
opened, we also made use of an ^^Fe calibration source 
located on the inside of the door, which illuminated the 
entire CCD imaging area. 



3.2. Energy scale offsets and origins 

Before describing the post-launch improvements made in 
our CCD response model, we address an important issue 
related to energy scale offsets, since they can lead to mis- 
interpretation of the data and hence strongly affect the 
modelling of the response, especially around the instru- 
mental edges. Below, we discuss four different causes re- 
sulting in energy scale offsets and describe the solutions 
which were found to correct the XRT energy scale. 

3.2.1. Evaluation of the bias level 

In addition to thermal noise, each pixel charge will carry 
with it a fixed zero-point offset DN value, a random read- 
out noise from the amplifier and the noise from camera 
electronics. The zero-point offset defines the bias level 
measured on the CCD. The amplifier readout noise and 
the electronic noise, however, affect the determination of 
the bias level, but we minimize those contributions by av- 
eraging multiple bias measurements. In any case, the noise 
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Table 1. Summary of the in-flight calibration targets used since the launch on 2004 November 20. 



Object^ 


Type 


Mode 


Purpose 


Exposure (ks) 


RXJ 1856.5-3754"'"'" 


Neutron star 


PG/WT 


Low energy response 


81/70 (31/40)^ 


PKS 0745-19 


Gluster of 
galaxies 


PG 


Effective area 


63 (20) 


2E 0102-723"^'"'" 


Line-rich SNR 


PG 


Gain, energy resolution and 


119 (28) 






WT 


shoulder 


76 (26) 


Gas A"'"'" 


Line-rich SNR 


PG 


Energy scale offset, gain, shoulder. 


235 (67) 






WT 


GTI and energy resolution 


62 (39) 


3G 273*/"^'^'"'" 


Quasar 


WT 


Effective area 


18 


PSR 0540-69"^ 


Pulsar 


PG/WT 


Effective area 


56/26 (25 in PG) 


PKS 2155-3047x 


Blazar 


PG/WT 


Effective area 


14/13 (5/9) 


NGG 7172*/" 


Seyfert 2 


PG 


Redistribution 


15 


G21.5" 


Featureless SNR 


PG 


High-energy shelf 


43 (78 in PG & 77 in WT) 


Mkn 421* 


Blazar 


WT 


Effective area 


33 (10) 


Grab"'"'"'"'^ 


Pulsar/SNR 


WT 


Effective area 


46 (5) 



The letter indices x, c, s, r, i correspond to sources observed by other X-/ Gamma- ray instruments for calibration purposes: 
XMM-Newton; Chandra; Suzaku; RXTE; INTEGRAL, respectively 

^ The numbers in parentheses correspond to the exposure time collected from September 2007 to the end of December 2007 



and dedicated to the re-calibration of the XRT after the substrate voltage change (see Section 4.3), while the other numbers 
correspond to the calibration data collected prior to the substrate voltage change on 2007 August 30. 

* This symbol corresponds to the sources for which simultaneous calibration observations were performed with the XMM- 
Newton EPIG cameras. 



Table 2. Summary of the releases of the XRT post-launch gain files. 



Release 
number^ 


Release 
date 


Main improvements in the gain 


Text section 
number 


05 


2005 April 5 


Introduction of a temperature-dependent gain file 


3.2.3 


06 


2005 October 28 


Update of the gain coefficients over time 




07* 


2007 July 9 


Update of the GTI values over time in WT & PG mode 

Update of the gain coefficients in WT & PG modes after the substrate 
voltage (Vss) change from OV to 6 V on 2007 August 30 


3.2.2 
4.3 


08 


2008 June 25 


Introduction of an offset of C3 = 17.6 eV to restore the WT energy scale 


3.2.4 



' The v003-v004 gain files correspond to the pre-launch gain files. 

The gain file has now added a suffix 'sO' or 's6' to indicate the appropriate substrate voltage (Vss)- The 'sO' files are valid 
from the start of the mission until when the substrate voltage was raised to Vss = 6 V on 2007 August 30. 



contributions should be much less than the measured off- 
sets. In orbit, bias frames in PC mode and bias rows in 
WT mode are taken during each slew of the spacecraft to 
a new target before the beginning of an observation. For 
PC mode, the bias is computed on a pixel- by-pixel ba- 
sis. Five 600 x 600 pixel bias frames are taken, averaging 
each into the existing bias map using a "running-mean" 
algorithm that sets each pixel (x, y) to a new mean value 
computed as follows: 

(N — 1) Pix fx v) 
MeanPixi+i(x,?/) ^ ^ ^ MeanPixi(x, y) + 

where TV = 3 is the running mean length and i = 1, 5 is 
the frame count. For WT mode, the bias row is a vector of 
bias values, one for each column in the WT mode window 



(i.e. 200 columns). The bias row vector is computed over 
one image of 600 rows. A "running- mean" approximation 
is computed for each pixel x in each row according to the 
following equation: 

MeanPixi+i(x) = ^ V eanPix^(x) + 

with = 3 and i corresponding to the row. The bias in- 
formation is then used to correct on-board the subsequent 
frames of data during the XRT observations. 

The bias level is mode dependent, and we have also 
seen observational evidence that the bias level can signif- 
icantly vary even during a single snapshot on a celestial 
target (i.e. over ~ 20 minutes). Variations of the measured 
bias level can be due to changes in the CCD tempera- 
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XMM: black - N, red - SE; XRT (bld26): green - N, blue - SE; 
XRT (bld27): cyan - N, magenta - SE 
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Fig. 9. (Top panel) Detector coordinate image showing 
the SNR Gas A when scattered optical light from the sunlit 
Earth is present on the detector (visible on the left side of 
the detector). (Bottom panel) PG mode light-curve of the 
SNR Gas A in the 0.3-10 keV energy range for different 
snapshots and grade selections: 0, 1-4, 5-12 and 13-31. The 
light-curve was obtained before bias correction. The drop 
of the count rate in the PG grade light-curve is due to the 
presence of scattered optical light from the sunlit Earth, 
which induces a grade migration from good grades (0-12) 
to higher grades. 



ture following some sky pointing directions. Shifts in the 
GGD bias level are also caused by roughly sinusoidal tem- 
perature variations of ±3°G during each 95 minute orbit. 
The measured bias level can also be contaminated by scat- 
tered optical light from the sunlit Earth as shown in the 
top panel in Fig. |9] (see Beardmore et al. 2007 for more de- 
tails). Scattered optical light from the sunlit Earth mainly 
affects PG mode because of its much longer frame accu- 
mulation time than WT mode, and seems to occur when 
the instrument is pointed near the sunlit Earth horizon. 
Ghanges in the bias level in both WT and PG modes result 
in energy scale offsets (see Fig. [lO|. 

In order to correct the bias level on the ground, and 
hence restore the energy scale (see Fig. 10), the task 
XRTWTCORR, for WT mode, was enabled in the xrtdas 
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Fig. 10. (Top panel) Gomparison of the event energy spec- 
tra of the Si and S lines in the North (N) and South 
East (SE) knots of the SNR Gas A as observed by the 
XMM-Newton MOS cameras (N: black; SE: red) and the 
Swift-XRT (N: green, cyan; SE: blue, magenta). The green 
and blue crosses correspond to XRT/PG grade data for 
which the bias was contaminated by optical light from the 
sunlit Earth. The data were processed with the v2.6 XRT 
software which does not correct for bias contamination. In 
this case, an energy scale offset is observed when compared 
with the XMM-Newton MOS spectra. The cyan and ma- 
genta crosses correspond to the same data processed with 
the XRT software including the task xrtpcbias, which 
corrects the data. (Bottom panel) Swift-XKT WT grade 
0-2 spectrum of RS Ophiuchi: (black) the data not cor- 
rected for the bias contamination and {red and blue) the 
data corrected using the modified task xrtwtcorr. 



software package (Swift software v. 2. 6). For PG mode, 
a new task xrtpcbias was developed and released with 
Swift software v. 2. 7. 

In addition to producing energy scale offsets, in some 
extreme cases scattered optical light in PG mode can 
induce a significant grade migration from good grades 
(grades 0-12) to higher rejected grades, resulting in an 
apparent loss of counts in the light-curve (see the bottom 
panel in Fig. l9|. In these extreme cases, the xrtpcbias 
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Fig. 11. (Top panel) WT grade 0-2 spectrum of the SNR 
Cas A in the energy band of the Si Ka and S Ka lines at 
different epochs since launch showing the degradation of 
the gain due to the CTI increase. The WT data were pro- 
cessed using the option biasdiff in the task xrtwtcorr 
correcting the bias level if necessary (see Section 3.2.1). 
The energy centroids of the Si Ka and S Ka lines show an 
energy shift of — 50eV from February 2005 to June 2007 
(see Section |4.2| for details about the line broadening). 
(Bottom panel) As above, but with the gain corrected for 
the effects of CTI. 



task will not be able completely to correct the bias level 
and a residual energy scale offset will be still present in 
the spectra. 

3.2.2. Increase of charge transfer inefficiency 

CCD detectors provide good X-ray imaging and spectro- 
scopic performance. However, the increase of charge trans- 
fer inefficiency (CTI) over time due to radiation damage is 
a fundamental limitation of CCD technology. Due to the 
effects of CTI, charge is lost during the readout process so 
that the remaining charge Q reaching the output amplifier 
after N, 

in the parallel direction is 



transfers in the serial direction and Np transfers 



Q = Qo(l-CTI,)^»(l-CTIj,) 



(9) 



where Qo is the initial charge and CTI^^^ are the serial and 
parallel CTI values, respectively. The CTI increase in the 




Fig. 12. (Top panel) Evolution of the gain Co at a CCD 
temperature of — 65°C. Note the jump in gain around 
MET (Mission Elapsed Time) 2.1 x 10^ s due to the 
substrate voltage change from Vss = OV to 6V (see 
Section [43] ) . (Bottom panel) Evolution of the serial (the 
blue line) and parallel (the green line) CTI values over 
time. Note the increase in CTI around MET 1.6 x 10^ s 
{i.e. January 2006). 



imaging area, the store-frame area and the serial register 
can then result in an energy scale offset if the data are not 



corrected. As an example. Fig. pT] shows the WT grade 0- 
2 Cas A spectrum in the Si Ka and S Ka line region at 
three different epochs when the data are not corrected for 
the effects of CTI. An offset of — 50eV is clearly visible 
between February 2005 and June 2007. 

The XRT gain file takes into account the effect of CTI 
as follows: 



PI = 



PHA X (Co + Ci X + C2 y) + C3 
G 



(10) 



where x and y correspond to the position of the event 
in the CCD imaging area. G = 10 is the global PHA to 
PI gain factor. Cq is the multiplicative DN to PHA gain 
factor, while Ci^2 = Cq x CTI^^^ are the multiplicative se- 
rial and parallel CTI correction factors, respectively (see 
Pagani et al. 2008 for more details). C3 is an additive off- 
set correction factor which was set to in both PC and 
WT modes before launch. Co, 1,2 are functions of the CCD 
temperature (see Section 3.2.3) and time. The CTI^^^ val- 



ues used to compute the Ci,2 coefficients in the gain file 
are the same for all CCD temperatures (see Section 3.2.3). 

Since launch, the evolution of serial and parallel CTI 
and the gain over time have been monitored using the four 
corner calibration sources (see Fig.|8|. The gain Co is given 
by the gain of the bottom left corner source CS3 closest 
to the output amplifier, as this does not suffer from CTI 
loss in the imaging area. Co thus includes both the output 
amplifier gain and the degradation in the gain caused by 
charge loss during the frame-store transfer. CTI5 is mea- 
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sured using corner source pair CS3 and CS4, while CTI^ 



measured using corner source pair CSi and CS3. Figure [12 
shows the evolution of the gain Co (from 2.529 just after 
launch to 2.625 in June 2007) and the CTlg^p values over 
time. Figure [12] shows a jump in gain Co around MET 
(Mission Elapsed Time) 2.1 x 10^ s due to the substrate 
voltage change from Vgs = OV to 6V (see Section 4.3). 



There also is a jump in CTI around MET 1.6 x 10^ s; the 
origin of this jump is however unknown. The correction of 
the data from the effects of CTI is illustrated in the bot- 



tom panel of Fig. 11 These CTI values were implemented 
in the PC and WT gain files from version 007 and onwards 
(see Table |2|. 



3.2.3. The effect of temperature 
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The XRT was designed to cool the CCD to a nominal op- 
erating temperature of -100°C using a thermo-electronic 
cooler (TEC). However, the TEC power supply system 
apparently failed shortly after launch, and the XRT has 
to rely on passive cooling via a heat pipe and radiator, 
in combination with enhanced management of the space- 
craft orientation to reduce the radiator view of the sunlit 
Earth. The XRT is now operated with CCD temperatures 
varying between -75 to — 52°C (see Kennea et al. 2007 for 
more details), though new GRBs occasionally cause it to 
exceed -52°C. 

Since the gain of the CCD output FET (Field Effect 
Transistor) is sensitive to the temperature, the reduced 
temperature stability required the introduction of a tem- 
perature dependency in the gain expression in order to 
restore the XRT energy scale. In addition to the Ka and 
KjS fluorescent lines of the corner source data, we used the 
supernova remnant Cas A, because its spectrum shows 
intense silicon (1.86 keV) and sulphur (2.45 keV) lines 
as well as an iron (6.6 keV) line. The gain coefficient 
Co in Eq. [lO] was computed at two CCD temperatures 
other than the — 100°C used during the ground calibra- 
tion: — 65°C and — 48°C using in orbit Cas A data. The 
ground data processing software linearly interpolates be- 
tween tables of gain coefficients for these three different 
temperatures to determine the correct value to use for a 
given observation!^ These temperature-dependent values 
of Co were used to compute the coefficients Ci,2 as well 
(see Section 3.2. 2[ ). The linear slope characterizing the gain 
change as a function of temperature from — 65°C to — 48°C 
is 0.00117 eV DN-^ °C-^ This corresponds to a - 7 eV 
variation in the energy scale at 1.5 keV for a 10° C tem- 
perature variation. This temperature dependency on the 
gain was included in the v005 release of the PC and WT 
gam files (see Table [2]). 
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Fig. 13. Spectra of the entire SNR Cas A in WT mode 
at different epochs: (black) 2005-02-16; (red) 2006-04-20; 
(green) 2007-06-06. The blue line indicates the expected 
energy centroid of the Si and S Ka lines. (Top panel) Data 
processed using the v007 WT gain file. (Bottom panel) 
Data processed using the v008 WT gain file including an 
offset of 17.6 eV. 



3.2.4. Comparison of the WT and PC energy scale 

We found that even after correcting the data for the var- 
ious sources of energy scale offsets discussed in previous 
sections, there was still a slight systematic mismatch be- 
tween the WT and PC energy scales (see the top panel in 



Fig. 13). We estimated the systematic offset to be 17.6 eV, 
based on the comparison of the line centroids between WT 
and PC Cas A spectra, as well as the use of the XMM- 
RGS model to fit the XRT/WT spectra of the SNR E0102- 
723 using the gain command within Xspec (Arnaud 1996). 
This energy shift appears to be independent of time, tem- 
perature and energy. The setting of the offset term C3 in 

Eq. [To] at 17.6 eV in the WT gain file was shown signif- 

^ http://swift.gsfc.nasa.gov/docs/heasarc/caldb/swift/docs/xrij?^antly to improve the WT energy scale (see the bottom 
SWIFT-XRT-CALDB-04_v2.pdf panel in Fig.fT3|. This new v008 WT gain file was released 
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on 2008 June 25 (see Table [2|. The origin of the energy 
offset in WT mode is unclear and still under investigation. 

3.3. Improvements to the response model 

We discuss in this section the improvements made to the 
CCD response model based on experiences from in- flight 
calibration. Table |3] summarises the improvements made 
to the response files and the corresponding release number. 

3.3.1. The low-energy response 

Shortly after launch, it became apparent that the low 
energy {E < 0.5 keV) response model could be im- 
proved. Figure [14] shows that the PC grade spec- 
trum of the soft neutron star RX J1856. 4-3754 (the 
black curve) is poorly fitted when the v007 RMF is 

used. A CONST*WABS*(BBODYRAD+BBODYRAD) model with 

the spectral parameters fixed to the values given in 
Beuermann et al. (2006) {i.e. Nh = 1.1 x 10^° cm-^ 
kTi = 62.8 eV and kT2 = 32.3 eV) was used for the 
fit. The lowest temperature black-body component has 
a minor impact in the XRT energy range. It was intro- 
duced by Beuermann et al. (2006) to fit the EUVE data 
as well as the Chandra data. Below 0.3 keV, the modelled 
line profile and its energy centroid are strongly dependent 
on the estimate of the charge losses at the CCD surface. 
As explained in Section 2.1.3, these charge losses are a 
function of energy and location of interaction of incoming 
photons. To model these charge losses better, we used the 
semi-empirical formalism described in Popp et al. (2000) 
to describe the loss function f{AEi=i^i2, z) in our CCD 
response model as follows: 

0, Z <0 OT z > AzccD 

/(Ai?i=i,i2,^) = { fo + A {fY, 0<z<l 



_ (z-l) 

Be - , I < z < AzccD 



where /o, c, / and r are free parameters which were es- 
timated initially by fitting the previous loss functions in 
each energy band AEi=i^i2 as defined in Section 2.1.3, 
and then by using an iterative process of RMF generation 
and fits of spectra of soft neutron stars and ground cal- 
ibration data. AzccD is the CCD thickness {Azccd = 
280 /im). The parameters A and B were derived using 
/(AEi=i,i2,z = I) and /'(AEi=i,i2,^ = I) = f U=;. This 
new description gives better results as shown in Fig. 14 
(the red curve). This description was included in the v008 
release of the PC and WT RMFs. 

3.3.2. The nnodelling of the shelf fronn photons above 
- 2 keV 

Before the release of the v009 RMFs, the spectral fits of 
heavily absorbed sources (with Nh typically larger than 
10^^ cm~^) in either PC or WT mode showed an underes- 
timation of the modelled redistributed counts at low en- 
ergies corresponding to the shelf (see the red curves in 




0.5 

channel energy (keV) 



Fig. 14. Comparison of the low-energy response be- 
tween vOOr (black) and v008 (red) RMFs using 
the spectrum of the soft neutron star RX J1856.4- 
3754 in PC mode for grade events. We used a 

CONST*WABS*(BBODYRAD+BBODYRAD) model with the 

spectral parameters fixed to the values given in 
Beuermann et al. (2006). 



Fig. 15). As discussed in Section 2.1.3, the physical ori- 



gin of the shelf, and hence its modelling, is complex. To 
improve the CCD response model, we decided to rescale 
the shelf for incident photons above 2 keV, since this is 
a straightforward approach. The result of this rescaling 
significantly improves the residuals at low energy when 
fitting spectra of heavily absorbed sources in both modes 



(see the black curves in Fig. 15 ). While the rescaling of the 
PC shelf did not change the QE at high energy because the 
shelf is at least two orders of magnitude below the photo- 
peak, the rescaling of the shelf in WT mode increased the 
QE by ~ 10% at 6 keV. This increase of the high-energy 
QE in WT mode also gives a better agreement with expec- 
tation of the grade ratios between the PC and WT modes 



(see Section 3.4). This new modelling was included in the 
release of v009 RMFs. 



3.3.3. Origin of the low-energy shoulder fronn 
high-energy photons 

For low-energy photons (below about 2 keV), the low- 
energy shoulder of the line profile can be modelled suc- 
cessfully by charge losses at the interface between silicon 
and oxide layer. However, the same process cannot com- 
pletely explain the low-energy shoulder seen in the line 
profile of high-energy photons. Hence, the pre-launch PC 
and WT RMFs (v006) were not able fully to model the 
line profile for high-energy photons (see the bottom panel 
in Fig.[7l). 

Tests performed after launch showed that an artifi- 
cial increase of the split threshold in the CCD response 
model results in an increase of sub-threshold losses, and 
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Table 3. Summary of the releases of the XRT post-launch RMFs. 



Release number^ 


Release date 


Main improvements in the RMFs 


Text section number 


07 


2005 April 5 


Line shoulder: ad-hoc increase of threshold to favour 
sub-threshold losses in PC & WT mode 


3.3.3 


08 


2006 April 24 


Low-energy response: change in the description the 
surface loss function in PC & WT mode 


3.3.1 


09^ 


2007 May 31 


Line shoulder: New description of the charge cloud in 
the field-free region for high-energy photons (E > 2 keV) 
in PC mode 

Shelf: Rescahng of the shelf in both PC k WT mode 


3.3.3 
3.3.2 


11* 


2008 June 25 


Loss function: New description of the surface loss function 
in the 1-2 keV energy range in WT mode 


3.3.4 



' The v006 response files correspond to the pre-launch response files. 

* The v009 response files were renamed as vOlO response files when the substrate voltage was raised from to 6 V on 2007 
August 30 (see Section 4.3) following an update of the ground software. 

* |http: / /swift .gsfc.nasa.gov /docs /heasarc /caldb /swift /docs /xrt /S WIFT-XRT-C ALDB-09- vl 1 .pdf (Godet et al. 2008) 
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Fig. 15. Best fits of absorbed sources {Nh > 5 x 10^^ cm ^) in PC and WT modes showing the improvement in 
the modelling of the shelf with the v009 RMFs. Top panel: (left) PC grade 0-12 spectrum of NGC 7172 using the 
v009 (black) and previous vOOS (red) response files. The model was a wabs*powerlaw model; (right) the same for 
the PC grade 0-12 spectrum of the SNR G21.5 using a tbabs*powerlaw. The spectrum was extracted using a 40" 
circle region centred on the core of the remnant. Bottom panel: (left) WT grade 0-2 spectrum of the X-ray binary 4U 
1608-52 using the v009 (black) and previous v008 (red) response files. The model was a wabs*(powerlaw+diskbb) 
model; (right) the same for the WT grade 0-2 spectrum of the microquasar GRO J1655-40 in a low/hard state using 
a WABS*(POWERLAW+GAUSS) model. 
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Fig. 16. Plot of the PC quantum efficiency using v007 
RMF for: grade 0-12 (red) and grade (black). 

hence a better modelling of the shoulder. This empirical 
technique was implemented in the v007 RMFs. Even if 
the global result was to improve significantly the residuals 
around the shoulder (see Fig. 6a in Osborne et al. 2005), 
this technique was not entirely satisfactory since residuals 
were still present and the physical origin was unidentified. 
The other drawback was that the QE curve in PC mode 
showed discontinuities (smaller than 5%) due to the use of 
different values of the threshold depending on the energy 



range (see Fig. 16), although these small discontinuities 
had no noticeable impact on the spectral fitting. 




channel energy (keV) 



Fig. 17. PC grade 0-12 spectrum of the ^^Fe source on 
the XRT camera door fitted using the v009 RMF. The low 
energy shoulders of the lines are now correctly modelled 
using the Pavlov & Nousek formalism (compare to Fig. [t] 
lower panel). 

Pavlov & Nousek (1999) have stressed that when a 
photon interacts in the field-free region, the resulting 
charge cloud no longer has the profile of a 2-D Gaussian, 



because it experiences free diffusion in this region be- 
fore being distorted when penetrating the depletion region 
where the charge is collected. The authors gave an ana- 
lytical formalism to describe the resulting shape of the 
charge cloud as a function of the depth of interaction. 
The resulting profile shows more extended wings which 
naturally favour the increase of sub-threshold losses for a 
given threshold value. The implementation of the Pavlov 
& Nousek formalism in our CCD response model gives 
good results as shown in Fig. [iTj This new modelling was 
included in the release of v009 RMFs. 

3.3.4. Feature around 0.9-1.0 keV in WT nnode 

Fits of high statistical quality WT spectra have revealed 
10% systematic residuals around 0.9-1.0 keV, as shown in 
Fig. [T|| An ad-hoc dip was added to the WT v008-v010 



ARFs around 0.9 keV in order to suppress these residu- 
als. However, this approach was not entirely satisfactory. 
These systematic residuals are due to a RMF redistribu- 
tion issue. Indeed, as explained in Section [2.1.3| the sur- 
face charge losses are modelled by a loss function / so that 
its parameters are different in different energy segments. 
1 keV corresponds to a boundary between two energy seg- 
ments for which there is a slight discontinuity in the de- 
scription of the loss function. We modified the parameters 
of the loss function in the 1-2 keV range in order slightly 
to change the monochromatic line input profile, as shown 
in Fig. [19] for 1.2 keV photons. New WT grade and 0-2 
RMFs {i.e. vOll) were computed using this new descrip- 
tion. These vOll RMFs give very good performance as 



shown by the red residuals in Fig. 18 (also see Section 4.1 ). 



3.4. Investigation of the QE shape from in-flight 
calibration 

3.4.1. Overview of the QE calibration 

The pre-launch QE was calibrated using measurements 
across the XRT energy range (0.3-10 keV) made at the 
Leicester calibration facility. However, some of these mea- 
surements were corrupted during the data taking process 
indicating that some of the low energy QE values were not 
valid. 

It is not straightforward to obtain a direct measure of 
the QE in orbit, since the overall instrument spectral re- 
sponse depends on the mirror, filter transmission and CCD 
responses. Early fitting results suggested that the effective 
area needed to be rescaled between modes and grade se- 
lection due to QE uncertainties. It was also necessary to 
introduce corrections around the silicon and oxygen edges 
to obtain flat residuals when fitting spectra in both PC 
and WT mode. As discussed in Section [3]2j the correction 
around the oxygen edge was eventually found to be needed 
due to energy offsets in both modes. 

We discuss below possible explanations accounting for 
the changes implemented in the ARFs and how they can 
be transfered to the CCD response model. 
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Fig. 18. WT grade 0-2 spectrum of Mkn421 from June 
2006 data fitted using the vOlO (black) and vOll (red) 
RMFs. The residuals around 0.9-1.0 keV (black points) 
are due to a RMF redistribution issue in v007-v010 WT 
RMFs. The use of the new vOll WT RMFs including a 
slight change in the loss function between 1 and 2 keV 
gives flatter residuals in this region (red points). The 
spectral model is an absorbed bending power-law with 
Nh = 1.6 X 10^0 cm-2. 




channel energy (keV) 

Fig. 19. Gaussian models with a value of a = 0.1 keV 
and an energy centroid of 1.2 keV folded through the RMF 
kernel: (black) CALDB WT grade RMF vOlO; (red) New 
WT grade RMF vOll. Note the change in the line profile 
due to a minor tweak of the loss function between 1 keV 
and 2 keV in order to remove the systematics observed 
around 0.9-1.0 keV. 



3.4.2. The silicon Ka edge and the low-energy QE 
shape 

Since the detector consists of silicon bulk, and the effective 
area has a maximum in the 1.5-2.2 keV energy band, it 
is essential to take special care when modelling the Si Ka 
edge (1.839 keV). 



Fits of several bright sources in both WT and PC 
modes using pre-launch RMFs demonstrated that the 
modelling of the silicon edge was not completely correct, 
as shown by the black residuals below 2 keV in Fig. [20] for 
WT mode. The residuals above 2.2 keV are due to an im- 
proper modelling of the Au M-shell edge region (2.0-3.5 
keV) in the theoretical ARF. This ARF is made of the 
filter transmission response swxftrans200ioioivoo5.fits 
and the effective area (v004), which is computed by a ray- 
tracing code (Cusumano et al. 2005). 




channel energy (keV) 

Fig. 20. Best fits of the WT grade spectrum of Mkn421 
using the WT grade vOll RMF together with the the- 
oretical ARF (black) and the new ARF vOll (red) which 
includes the correction of the CCD QE just above the Si 
Ka edge (see text in Section 3.4.2). The residuals seen in 
black above 2.2 keV are due to the improper modelling 
of the Au M-shell (2.0-3.5 keV) edges in the theoretical 
ARF. A correction in the Au edge region (2.2-4 keV) was 
introduced in the ARFs after launch in order to flatten 
the residuals above 2.2 keV (red residuals). 



Initially, a correction was implemented in the v007- 
vOlO PC and WT ARFs to flatten the residuals in the 
1.8-2.2 keV energy range. However, that correction suffers 
from the effects of energy scale offsets which could not 
be entirely corrected at those times. Since then, we have 
been able to correct the XRT data from these effects (see 
Section 3.2). This enabled us to establish that the problem 
was related to an over-prediction of the QE in our CCD 
response just above the Si Ka edge (1.84-2.2 keV) rather 
than a redistribution issue. A straightforward solution to 
correct the QE above the Si edge was to implement a cor- 
rection based on the XMM-MOS QE curve, which shows a 
deeper profile at these energies. Unlike the Swift-XBT^ the 
XMM/MOS QE from the XMM-SAS 7.1.2 was calibrated 
with a synchrotron beam in Saclay, allowing for accurate 
measurements around the Si edge. This correction, along 
with a correction of the Au edge region (2.2-4 keV) in the 
ARFs produced after launch, gives good results in WT 
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mode (see the red curve in Fig. 20). A similar correction 
was implemented in the new vOll PC ARFs. 

Keay et al. (1995) stressed that the QE shape just 
above the Si edge was strongly dependent on the electrode 
structure and composition of their JET-X CCD. They 
showed that the silicon contributes 65% of the QE shape 
above the edge, while the oxide and nitride of silicon con- 
tribute at 28% and 7%, respectively. As shown in Fig. [s) 
the description of the electrode structure in the CCD pixel 
geometry used to generate the XRT RMFs is relatively 
simple. In reality, there is an overlap of three different 
electrodes. An explanation for the QE over-prediction just 
above the silicon edge could be that the thicknesses of the 
different materials in the CCD geometry need to be more 
accurately specified. However, variation in thickness of the 
Si and Si02 layers in the electrode and finger structures 
did not allow us to converge towards adequate thickness 
values for these layers to fully model the QE just above 
the silicon edge. 

Nevertheless, from in-flight calibration we noticed that 
the total effective area at low energy needs to be decreased 
in both WT and PC modes to obtain correct flux lev- 
els and spectral parameters. Again, the electrode thick- 
nesses appeared to be under-specified, since more low- 
energy photons will be stopped in the finger and electrode 
structures. 

Another possible explanation accounting for the QE 
reduction at low-energy is that the etching of the elec- 
trode and finger structure did not result in the designed 
electrode widths, and some residual material was left as 
shown in Fig. [3] The residual material acts as an extra 
absorption layer for low-energy photons and so reduces 
the QE at low energy. As a test we reduced the open elec- 
trode area by slightly increasing the channel stop width 
(by 1 /im) and the finger as shown in Fig. |2l] (top panel), 
corresponding to an open area reduction of 10.6%. The 
bottom panel in Fig. |2l] shows a comparison of the WT 
grade 0-2 QE curve when using the new geometries (in red 
and green) and that from the WT grade 0-2 RMF vOlO 
(in black). A less than 10% decrease of the QE is obtained 
at the energy of the oxygen edge (0.543 keV). This is con- 
sistent with the ad-hoc reduction implemented in the WT 
v007-v010 ARFs. 

New WT ARFs (vOll) were created by rescaling the 
theoretical ARF below 2 keV using the green QE curve in 



the lower panel in Fig. [21] and by introducing the CCD QE 
correction just above the Si Ka edge as discussed above, 
while the Au M-shell region was corrected using celes- 
tial continuum sources. These new WT grade and 0-2 
ARFs have similar shapes except in the 1.5-1.84 keV en- 
ergy range where the WT grade 0-2 ARF is rescaled down 
by less than 2% with respect to the WT grade ARF in 
order to flatten residuals below 2% (see Fig. 3 in Godet 
et al. 2008). Figs. [20] (red curve) & [23] show WT spectra 
of Mkn421 fitted using an absorbed bending power-law 
and the new vOll response files. This plot shows the very 
good performance of the vOll WT RMFs/ ARFs, since 
the residuals across the 0.3-10 keV energy range are bet- 
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Fig. 21. (Top panel) Drawing showing the changes in the 
geometry of the finger and channel stop, corresponding 
to a reduction of the open area by 10.6%. (Bottom panel) 
Plot of the low-energy quantum efficiency for different con- 
figurations of the geometry of the CCD pixel: (black) using 
the default CCD geometry shown in Fig. [s] (red) with an 
increase of the width of each channel stop by 1/im; (green) 
with an increase of the width of each channel stop by 1/im 
and the finger structure as shown in the top panel. 

ter than 3% and the systematic errors are less than 2% 



(also see Section 4.1). 



3.4.3. Discrepancy in the effective area between PC 
and WT nnodes 

We expect that differences in the effective area between 
PC and WT modes would be due to the differences in 
the event grading in these two modes. However, it was 
necessary to scale down the effective area in PC mode (for 
both grade and 0-12) by a factor larger than expected 
from the simple differences in the event grading, in order 
to reproduce the correct flux levels for the same calibration 
targets (see the top panel in Fig. 



22 



see Section 3.2 in 
Godet et al. 2008). Indeed, while the shape of the PC 
grade 0-12 effective area is similar to that of the WT grade 
0-2 effective area above the Si Ka edge (as expected from 
the event grading), the WT grade effective area is larger 
than the PC grade 0-12 effective area below the Si Ka 
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edge. The event grading in PC and WT mode implies 
the opposite, since WT grade events only contain PC 
grade 0, 1 & 3 events. The reason why the PC grade 0-12 
effective area is less than the WT grade effective area 
below ~ 1.84 keV remains unclear. 



WT grade / 
WT gi ades 0-2 / 
PC gi ade ' 
PC gi-ades 0-12 




Chaimel energy (keV) 



Fig. 22. (Top panel) Comparison of the XRT total ef- 
fective area (EA) in PC and WT mode using the vOll 
response files. The PC grade EA is very low when com- 
pared to the PC grade 0-12 and WT EA even at low en- 
ergy {E < 1 keV) where we expect to have a similar QE 
level whatever the mode and grade selection. Below ~ 1.84 
keV the PC grade 0-12 EA is less than the WT grade 
EA, whereas we expect the WT grade events to con- 
tain the PC grade 0, 1 and 3 events. (Bottom panel) PC 
(black: grades + 1 + 3) and WT (red: grade 0) spectrum 
of the whole remnant Cas A. The PC and WT data were 
processed using a central threshold of 80 DN and a split 
threshold of 80 DN. For both spectra, we used the same 
size of the extraction region. There is a clear discrepancy 
between the two spectra between 1 and 3 keV which can- 
not be due to differences in the response redistribution 
between PC and WT mode. 



There is also a difference in the effective area at low 
energy (below 0.5 keV) between the PC grade and WT 
grade effective areas. The fact that the WT grade 
events contain the PC grade 0, 1 & 3 events may account 
for the discrepancy. However, the lower panel in Fig. 22 
shows that there are some differences in the Cas A spectra 



between WT grade and PC grades + 1 + 3 in the 1-3 
keV band. To extract these spectra, we used 2005 data col- 
lected prior to the apparition of the bad columns, and the 
data were processed using the same central (80 DN) and 
split (80 DN) event thresholds. Below - 1 keV, the differ- 
ences observed between the two spectra may be partially 
due to differences in the kernel redistribution between the 
modes (see Fig. |6|. That suggests that the observed dif- 
ferences in effective area at low energy cannot be due to 
simple differences in the event grading. 

4. In flight spectroscopic performance and caveats 

We discuss in this section the in-flight spectroscopic per- 
formance of the vOll response files when fitting data col- 
lected at a substrate voltage of Vgs = OV. For all the 
fits shown below, the PSF and exposure map corrections 
were taken into account in the ARFs using the version 
5.5 of the task xrtmkarf of the xrtdas software package 
(Capalbi et al. 2005). We compare the fitting results with 
those found using other in-orbit X-ray instruments. We 
also discuss current caveats in the XRT spectral response: 
i) the line broadening due to the build-up of charge traps 
on the CCD over time and methods to handle it; ii) the 
impact of the permanent increase of the substrate voltage 
from V to 6 V on the performance of the vOll response 
files. 



4.1. Spectroscopic performance 

Figure |23] (panel a) shows the best fit of the high statis- 
tical quality WT spectra of the blazar Mkn421 observed 
in June 2006 using the vOll response files. Each spectrum 
contains more than 10^ counts. The residuals across the 
XRT energy range are below 3%. The fits of high statis- 
tical quality PC spectra show that residuals are less than 
5% across the XRT energy range (see Godet et al. 2008). 

The existence of several X-ray instruments (XMM- 
MOS, XMM-PN, Chandra, RXTE-FCA, Suzaku-XlS) 
covering an energy band similar to that of the XRT en- 
abled us to perform cross-calibration using different con- 
tinuum and line sources (see Fig. [23| panels b, d & e and 
Table |4j see Plucinsky et al. (2008) for a cross-calibration 
work on the SNR E0102-723). Figure |23] shows that the 
XRT spectral fit residuals are comparable with those of 
other instruments. The spectral parameters derived using 
XRT data are consistent with those derived using XMM- 
Newton, Chandra and RXTE dsid. (see Table|4|. The XRT 
fluxes from PC grade events seem to be slightly higher 
up to 7% when compared to those found from XRT PC 
grade 0-12 and other instruments, indicating that the PC 
grade effective area may need to be rescaled in the fu- 
ture. The WT grade and grade 0-2 fluxes derived from 
the fits of our WT calibration sources agree with each 
other within the error bars (see Table [4|. 

The spectral analysis of joint BAT-XRT GRB spectra 
during overlapping time intervals showed relatively good 
agreement between the two instruments, within 10% in 
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Fig. 23. Plots showing the best fits to different celestial sources when using the vOll PC and WT response files, (a) 
WT grade (black) and grade 0-2 (red) spectra of the blazar Mkn421 observed in June 2006. The spectral model is 
an absorbed bending power-law (i.e. wabs * exp[-a+/3 ln(E)]) with the column density fixed at Nh = 1.6 x 10^^ cm~^ 
(Kalberla et al. 2005). Each spectrum contains over 1.6 x 10^ counts. The spectra were heavily rebinned to better 
emphase the statistical residuals obtained for such high statistical quality data. The residuals are better than 3% over 
the XRT energy range, (b) Simultaneous XMM-Newton, Swift and RXTE observations of 3 C273 performed in July 
2005 showing the PN (black), MOSl (red), M0S2 (green), RXTE (blue) and XRT WT grade 0-2 (cyan) fit using a 
WABS*POWERLAW + 2 zBBODY model, (c) Simultaneous Swift-XRT and Swift-B AT observation of Sco X-1 performed 
in April 2008 showing the XRT WT grade 0-2 (black) and BAT (red) fit. The data were fitted using a const*wabs 
(bbodyrad + comptt) model with const fixed to 1 for the XRT data and left free when fitting the BAT data (const 
= 0.93 ± 0.04, the errors being given at 2.71 a). The XRT/WT spectrum was extracted using an annulus region with 
an inner radius of 25 pixels since the data were heavily piled-up due to the high count rate of the source (2020 counts 
s~^ when using a full circle extraction region), (d) Fits of PSR 0540-69 (the pulsar and the nebula): (blue) Swift-XBT 
PC grade 0-12; (green) XMM-PN pattern 0-4; (red) XMM-M0S2 pattern 0-12; (black) XMM-MOSl pattern 0-12. 
The model is a wabs(vnei + powerlaw) model, (e) Fits of the cluster PKS 0745-19: (red) Swift-XRT PC grade 0-12; 
(black) XMM-MOSl pattern 0-12. The model is a wabs*mekal model with the abundances left free. 
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most cases. Figure [23] (panel c) shows the joint fit of the 
BAT and XRT data of Sco X-1. The data were fitted using 
a coNST*WABS(BBODYRAD+coMPTT) model. The constant 
factor (const) was fixed to 1 for the XRT data and left as 
a free parameter for the BAT data (const = 0.93 ± 0.04, 
the errors being given at 2.71a). The residuals observed 
for the XRT spectrum in Fig. |23] (panel c) above 8 keV cor- 
respond to Ni and Au L-shell fluorescence. Discrepancies 
larger than 10% between the BAT and XRT flux normali- 
sation and/or large residuals at high-energy (above 5 keV) 
may mean that: i) the spectral model used is not suitable 
to fit the data; ii) the XRT data may be piled-up (pile-up is 
dependent on the spectral shape of the source. However, a 
rough limit for pile-up to become an issue is 100 counts s~^ 
in WT mode and 0.6 counts s~^ in PC mode for powerlaw- 
like spectra); iii) instrumental lines such as Ni Ka & K/3 
and Au L-shell fluorescence lines may be present above 8 
keV (see Moretti et al. 2008b). 

Our current understanding of the XRT response at 
Vss = OV implies a systematic error of less than 3% in 
both WT and and PC mode in the 0.3-10 keV energy 
band and better than 10% in absolute flux (see Godet et 
al. 2008). 

4.2. Line broadening 



In Section |3.2.2[ we discussed the effect of the increase 
in CTI (due to charge traps too shallow to be individu- 
ally identified) which can cause an energy scale change if 
not corrected. Besides this CTI, there is also a build-up 
of deeper charge traps due to high-energy proton and ra- 
diation damage to the CCD (the imaging area, the store 
frame area and the serial register), which dominate the 
line broadening observed in XRT data (see Fig. 11). The 
FWHM at 1.86 keV in the Gas A spectra has degraded 
from 105 eV in February 2005 to 131 eV in June 2007. 
Over the same interval, a 44% width increase was mea- 
sured at 5.9 keV from the on-board calibration sources. 
The most serious of these charge traps can cause a loss of 
up to 350 eV from the incident X-ray energy. 

Since launch, regular observations of the soft neutron 
star RX J1856. 4-3754 have enabled us to monitor any sig- 
nificant low energy redistribution change in the CCD re- 



sponse. The top panel in Fig. 24 shows the evolution of 
PC grade spectra of RX J1856. 4-3754 between February 
2005 and June 2007. Both spectra were obtained using 
time intervals when the source was not located on the bad 
columns. Fitting the spectra with the model described in 
Fig.[l4]gives a constant factor of CON ST = 0.97±0.03 for 
February 2005 data and CONST = 0.76 ± 0.02 for June 
2007 data. The use of a simple wabs*(bbodyrad) with the 
column density fixed at 1.1 x 10^^ cm~^ gives a black-body 
temperature of kT = 62.8±1.7 eV for both datasets; which 
is consistent with the canonical value found in the litera- 
ture. A similar evolution is also observed in WT mode with 
a constant factor of CONST = 0.98 ± 0.03 for February 
2005 data and CONST = 0.80 ± 0.03 for June 2007 data. 
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Fig. 24. Effects of the build-up of charge traps on the 
XRT spectral response for continuum sources. (Top panel) 
Evolution of the PC grade spectrum of the soft neu- 
tron star RX J1856. 4-3754 over time: (black) data from 
February 2005 and (red) data from June 2007. The de- 
crease in the flux normalisation between the two epochs 
by ~ 22% is due to the loss of events below the on- 



board central event threshold (see text in Section 4.2). 
(Bottom panel) WT grade (black) and 0-2 (red) spec- 
tra of Mkn421 observed in March 2007 fitted using the 
vOll WT RMFs/ARFs. The spectra were fitted using an 
absorbed bending power-law as defined in Fig. [23] The 
WT spectra contain more than 2 x 10^ counts. The resid- 
uals around the instrumental edges are due to the energy 



shifting effect of charge traps (see text in Section 4.2). 



All the errors quoted above are given at 2. 71 a. Continuum 
spectra of Mkn 421 from March 2007 WT data also showed 
spectral fit residuals around the instrumental edges, espe- 



cially the oxygen edge (see the bottom panel in Fig. 24), 
while no evidence for such residuals was observed in June 
2006 data. 

All these changes are likely to be due to the build-up 
of charge traps over time. Charge traps are thought to 
be due to faults in the silicon crystalline structure which 
trap a fraction of the charge passing through them during 
the readout process. Therefore, only the events occurring 
above a trap in a given column experience a charge loss. 
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The top panel in Fig. 25 gives an example of a charge 
trap formed just after launch in the column DETX 78. 
The operation of the CCD at relatively high temperatures 
(from -70°C to -50°C) may partially fill the shallowest 
traps since the level of the thermally induced dark current 
is higher than that expected at the designed operating 
temperature of — 100°C (see Section 4.3). 
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Fig. 25. (Top panel) PHA distribution of the Mn Ka and 
Kp line PC grade 0-12 events over the column DETX 78 of 
the CCD-22. The location of a charge trap is clearly visible 
around row 271. Only the events above the trap, i.e. above 
the row 271, are affected. (Bottom panel) WT grade 0-2 
spectra of the SNR Cas A: (black crosses) the spectrum 
from June 2007 data; (red crosses) the same when the data 
were corrected from the effect of deep charge traps by ap- 
plying energy offsets in different columns; (green crosses) 
spectrum from February 2005 data. The FWHM improves 
significantly after the correction of the data. 



Trap-induced line broadening is difficult to model since 
it may depend on the energy of the incident photons and 
the source intensity. Another uncertainty comes from the 
fact that the location and depth of the charge traps are 
unknown. Because all XRT spectral response files released 
so far have an incomplete modelling of the effect of CTI or 
traps, we urge caution in the interpretation of apparently 
sharp spectral features observed in XRT data from March 
2007 onwards using the current calibration files. 



Prospects to deal with the line broadening are under 
investigation: 

Characterization of traps with the largest depth - We 
are experimenting with a technique based on the char- 
acterization of the location and energy offset induced by 
charge traps with the largest depth. The bottom panel in 



Fig. [25] illustrates how powerful this technique can be in 
restoring the spectral resolution of the XRT; the Cas A Si 
line FWHM measured from corrected WT data from June 
2007 is 109 eV compared to 131 eV before correcting the 
data (see Godet et al. 2007b). 

Broadening of the RMF kernel - We are developing a 
new tool which allows us to broaden the RMF kernel by 
convolving each of the 2400 RMF spectra with a flux- 
conserving function / with a CTI-dependent FWHM. All 
our experiments so far focused on WT mode, but the prin- 
ciple should be similar in PC mode. The line broadening is 
dominated by the effect of the deepest charge traps, which 
shift a given spectrum to lower energy by some amount. 
It is also important to take into account the fact that the 
amount of charge lost depends on energy in a subtle and 
poorly understood way. 

Since the true energy dependency is still unknown, we 
divided the 0.3-10 keV energy range into three ranges: 0.3- 
2 keV, 2-5 keV and 5-10 keV. Then, we defined in each of 
the three ranges an ad-hoc shape of the function /; i.e. in 
each range we assumed that the shift induced by a trap 
is constant and independent of the photon energy. The 
function / is defined as the sum of two Gaussians, with 
their width, relative normalisation and separation kept as 
free parameters to be optimised using data from celestial 
sources, so that the FWHM of each RMF monochromatic 
line photo-peak of the vOll WT RMFs, once convolved 
by /, is multiplied by a factor B which tracks the CTI 
increase over time. 

To test this technique, we first calibrated the function 
/ using March 2007 data, and we computed experimen- 
tal WT RMFs with a broadened kernel. As illustrated in 
Fig. [26] (the top left panel), which shows the best fits of 
the March 2007 WT spectra of Mkn421 (with an average 
count rate of 45 count s~^) using these RMFs along with 
the vOll WT ARFs, the residuals around the instrumen- 
tal edges are flattened when compared to Fig. [24] (bottom 
panel). We also investigated whether variation of the CCD 
temperature could have a significant impact on the trap 
properties, and so the line broadening. To do so, we con- 
sidered four temperature ranges with a width of 5°C from 
-70° C to -52° C using Mkn421 data from March 2007 up 
to August 2007. We did not notice any significant change 
in the residuals around the instrumental edges when fit- 
ting the WT spectra using the experimental WT RMFs, 
even though the monitoring of the on-board corner sources 
has revealed a trend with temperature (a ~ 25% increase 
in the FWHM at 5.9 keV from -70°C to -50°C). 

We also tested these experimental WT RMFs on data 
from June 2007 up to September 2007 for other celestial 
sources with lower count rates, such as the SNR E0102- 



Godet et al.: Modelling of the Swift-XRT spectral response 



25 




0.1 



1 2 
channel energy (keV) 




+ 



0.5 

channel energy (keV) 




Energy (keV) 



Fig. 26. WT spectra of celestial sources using the vOll WT ARFs and/or the experimental broadened WT RMFs: 
(Top left panel) Fits of the Mkn421 data from March 2007 using the experimental broadened WT RMFs. The WT 
grade (black) and 0-2 (red) spectra contain more than 2 x 10^ counts and were fitted using the same model as given in 
Fig. [24] (bottom panel). (Top right panel) SNR 2E 0102-723 from June 2007 data (black: vOll WT grade 0-2 RMF; red: 
new broadened WT grade 0-2 RMF). The WT grade 0-2 spectrum contains more than 2.2 x 10^ counts. The spectral 
model is based on an XMM/RGS-PN model (see Plucinsky et al. 2008). The spectrum was extracted using data when 
the source was not located on the bad columns. We did not apply any PSF correction since the source is extended, 
and there is no tool in the ground software to generate extended ARFs. (Bottom panel) RX J1856. 4-3754 from June 
2007. The spectra contain more than 2.9 x 10^ counts. The spectral model is the same as given in Section 3.3.1 The 
value of the constant factor is C = 0.90 ± 0.02 for grade (black) and C = 0.94 ± 0.03 for grades 0-2 (red). All the 
data were processed with the new WT gain file swxwtgainso_200ioioiv008.fits (see Section 3.2.4). 



723 (1 count s~^) and the soft neutron star RXJ1856.4- 
3754 (0.3 count s~^), to check the RMF kernel broaden- 
ing and for any temporal evolution of the line broadening. 
Indeed, the higher the source count rate, the more likely 
the charge traps are to be filled; thus, it could be ex- 
pected that brighter sources show a less significant broad- 
ening. However, since the XRT operates at much higher 
temperatures than it was designed for, most of the small 
charge traps should be already filled. Figure [26| shows the 
very good performance of these experimental WT RMFs. 
Specifically we see no intensity effect for sources between 
0.3 count s~^ to 70 count s~^. Moreover, their use en- 
ables us to retrieve in the case of RX J1856. 4-3754 (June 
2007 data) a normalisation consistent with that found us- 
ing 2005-2006 PC & WT data (see Section [3XT| . This 
is because the new WT RMFs show a slight decrease of 
the QE at low energy to take into account the fraction of 



events lost below the on-board event threshold due to the 



effect of charge traps (see Fig. 27). 



4.3. Increase of the substrate voltage to 6V 



As discussed in Section |3.2.3[ the XRT operates at 



higher than expected temperatures, resulting in signif- 
icant thermally-induced noise appearing as low energy 
events. From experiments performed at the Leicester cal- 
ibration facility on the Swift flight spare CCD-22 devices, 
we have demonstrated that raising the substrate voltage 
to Vss = 6 V reduces the dark current, since the volume of 
silicon in which carriers are generated is reduced (Osborne 
et al. 2005, Godet et al. 2007b). Figure |28] clearly illus- 
trates the benefit of raising the substrate voltage, since, 
for the same level of dark current, it is now possible to op- 
erate the CCD at a 3-4° C warmer temperature, and hence 
to collect useful science data up to -50° C, before excessive 
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Table 6. Best fit parameters of the PC grade 0-12 spectra 
extracted using the 40 arc-second core of the SNR G21.5. 



Vss 


Nh 


r 


Normalisation^ 


(V) 


(xlO^^ cm-2) 




(xlO"^ 
ph keV"^ s-^) 









6 


q qO + 0.14 

0.00_Q 23 


1.92 ±0.07 
2.01 ±0.06 


-, 7-1+0.20* 
' -"--0.18 



^ The fits were performed using a TBABS*POWERLAW 
model with the abundance table given by Wilms et al. 
(2000). No PSF correction was applied since the source 
is extended. 

* All the errors quoted above are given at 2. 7 la. 



Fig. 27. WT grade 0-2 QE curves in the 0.2-1 keV energy 
band: (black) RMF vOlO, (red) RMF vOll and (green) 
broadened WT RMF vOll. The loss of QE at low energy 
in the broadened WT RMF mimics the loss of events below 
the on-board central energy threshold due to the effect of 
charge traps. 

hot pixels compromise the data. The only minor drawback 
of raising this voltage is that it induces a decrease of the 
depletion depth, and therefore of the QE at high energy 
and just below the Si edge (~ 1.5 — 1.84 keV), due to the 
lower transparency of the silicon in this energy range (see 
Fig. 4). The increase of the substrate voltage also results 
in a change in the gain Co, since the gain of the output 
FET is also modified in this new configuration. 



Numerical simulations performed using our CCD re- 
sponse model and laboratory measurements made on a 
spare detector have shown that the QE decrease should 
be small (less than 10% at 6keV; Osborne et al. 2005, 
Godet et al. 2007b). Short observations of two bright ce- 
lestial targets with Vss set to 6V - Gas A (2.2 ks in PC 
mode and 0.8 ks in WT mode) and the Crab (0.6 ks in WT 
mode) - were also performed to estimate the on-board ef- 
fective area change and to create new PC and WT gain 
files. The QE changes were measured to be ~ 7% at 6 
keV and less than 10% in the 1.5-1.84 keV range. These 
Vss = 6V gain files were released as version 007 prior to 
the permanent change (see Table [2|. The change in sub- 
strate voltage has made it necessary to release two sets 
of gain files, now distinguished by the characters 'sO' and 
's6' in their file names. 




Fig. 28. Evolution of the thermally induced dark current 
as a function of the CCD temperature: (blue points) before 
the substrate voltage change to = 6V; (red triangles) 
after the change. The level dark current (given in DN) 
is computed as the difference between the mean of the 
noise peak and the mean of the bias level in the data. A 
clear reduction of the dark current can be seen due to the 
reduction of the volume of silicon in which carriers are 
generated is reduced. 



Fig. 29. Comparison of the residuals for PC grade 0-12 
data of the SNR G21.5 collected at = (black) and 
6 V (red). The PC spectra were extracted using the 40 
arc-second core of the remnant. The spectra were fit using 
the vOll PC RMFs and ARFs and the v008 PC gain files. 
The PC spectra at Vss = and 6V contain 9.4 x 10^ 
and 1.5 X 10^ counts, respectively. The level of residuals is 
consistent between the two datasets as well as the spectral 
fitting parameters (see Table |6|. 
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Because the QE reduction was estimated to be rather 
small and the effect of the operational change on the spec- 
troscopic performance was demonstrated not to be signif- 



icant for most of the XRT observations (see Fig. 29 and 
Table |6|, the substrate voltage was permanently raised to 
6 V on-board on 2007 August 30. 

Even though the QE change is small, and really no- 
ticeable only in WT spectra with more than 10^ counts, 
the RMFs and ARFs need to be updated for both PC 
and WT modes. An intense phase of re-calibration of the 
instrument is ongoing, as indicated in Table 1. The full de- 
tails of the calibration of the Vss = 6V response files will 
be addressed in a forthcoming paper. By the time this pa- 
per is published, new response matrices will be available 
that take the QE change into effect properly. 

5. Conclusion 

We described in detail our Monte-Carlo simulation, com- 
puting the Swift-XRT PC and WT RMFs. The response 
model is mainly based on a physical description of the in- 
teraction of photons in the CCD. We showed how we used 
in-flight calibration to improve the XRT spectral response 
(the low-energy response, the line profile and the shelf) by 
implementing empirical corrections when it was not pos- 
sible to implement physical ones. All the changes allow us 
to describe the CCD response well. We showed that the 
vOll XRT response files, calibrated using data collected 
at Vss = V, give good performance on continuum and 
line sources in both PC and WT mode when compared 
to other X-ray instruments in the 0.3-10 keV energy band 
(the recommended bandpass) with a systematic error of 
less than 3% in both modes over 0.3-10 keV and better 
than 10% in absolute flux. 

We also showed that the XRT spectral response cal- 
ibration was complicated by various effects, such as en- 
ergy offsets, related to the way the CCD is operated in 
orbit. We described how these effects can be corrected 
in the ground processing software thanks to the task 
XRTWTCORR in WT mode and the task xrtpcbias in PC 
mode, and the release of temperature-dependent gain files. 
We stressed that the CCD started showing evidence of 
degradation of its spectroscopic performance (line broad- 
ening, small residuals around the instrumental edges and 
change in the low-energy response) when fitting spectra of 
celestial sources from the middle of 2007. This is due to the 
build-up of charge traps in the CCD, produced by high- 
energy proton and radiation damage. We are developing 
RMFs with a broadened kernel, enabling us to handle the 
line broadening over time. The preliminary results in WT 
mode are promising. Before the release of broadened ker- 
nel RMFs in either mode, we recommend that users should 
be cautious in any interpretation of the data. 

In order to improve further the XRT spectral perfor- 
mance, the substrate voltage was raised permanently to 
Vss = 6V on-board on 2007 August 30. This change re- 
sulted in a slight decrease of the QE at high energy {E > 3 
keV) and just below the Si edge (1.0 - 1.8 keV). This QE 



change is really only noticeable when fitting high statisti- 
cal quality spectra with more than 10^ counts. An intense 
phase of re-calibration of the instrument is on-going to 
update the spectral response files. In the meantime, we 
recommend the use of the vOll response files, and we ad- 
vice the user to be cautious in the spectral analysis of XRT 
data until the response files are updated. 
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